

%% specify files

matfiles = {'run1a_baseline_sr.mat','run1b_baseline_no_sr.mat',...
    'run5a_crra_sr.mat','run7a_Nstate_sr_fit_v3.mat','run6a_time_varying_s_sr.mat'};

for i = 1:5
    load(matfiles{i})
    stats_list{i} = stats;
end

%% term structure across models

for i = 1:numel(matfiles)
    load(matfiles{i})
    % nominal yields
    fprintf('\n\n%s, nominal:\n\n', matfiles{i})
    fprintf('%2.2f &', stats.nom_yld_mean)
    fprintf('& %2.2f', stats.nom_yld_std)
    % real yields
    fprintf('\n\n')
    fprintf('%s, real:\n\n', matfiles{i})
    fprintf('%2.2f &', stats.real_yld_mean)
    fprintf('& %2.2f', stats.real_yld_std)
end